Half-skyrmion picture of single hole doped high-T c cuprate 
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The single hole doped CuCh plane is studied by the 0(3) non-linear a model and its CP rep- 
resentation. It is argued that the spin configuration around a Zhang-Rice singlet is described by a 
half-skyrmion. The form of the excitation spectrum of the half-skyrmion is the same as that in the 
7r-flux phase. 
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In the phase diagram of high-temperature supercon- 
ductors, the most established phase is apparently the 
antifcrromagnetic long-range ordered phase, or the Neel 
ordered phase of undoped parent compounds. Due to a 
large charge-transfer gap [[J, a hole at each copper site in 
the CuC>2 plane is localized to form a spin S =1/2 mo- 
ment. There is an antiferromagnetic superexchange in- 
teraction J between the copper site spins and the system 
is described by the antiferromagnetic Heisenberg model 
on the square lattice, J¥afh = JJ^/i j\ Si • Sj. Antiferro- 
magnetic spin wave dispersions observed by neutron scat- 
tering experiments are in quite good agreement with the 
spin wave theory with the value of J determined by Ra- 
man scattering pj . High-temperature superconductivity 
occurs upon moderate hole doping on the Neel ordering 
phase. 

For understanding the mechanism of high-temperature 
superconductivity, it is necessary to describe the doped 
holes. However, for the purpose of figuring out the proper 
picture, considering a moderately doped system does not 
appear to be promising because such a system is quite 
complicated and indeed there is no consensus on the de- 
scription of the system. Contrastingly, a system with 
only one doped hole is much simpler and there seems 
to be much hope to understand the physics because it 
is closest to the well-established antiferromagnetic long- 
range ordered phase. Experimentally angle resolved pho- 
toemission spectroscopy (ARPES) 0, HQ on undoped 
compounds provide valuable information on the single 
hole doped systempj- 

As an effective theory of the antiferromagnetic Heisen- 
berg model Hafh, the 0(3) non-linear a model (NLerM) 
has been studied extensively ||. The model is derived 
from Hafh by applying Haldane's mappingQ: 



S = 



dr / d 2 r 



dn 



(1) 



where p® is the bare spin stiffness constant on the length 
scale of the lattice constant, which is taken to be unity, 
and c sw is the spin-wave velocity. The integral with re- 
spect to the imaginary time r is carried out up to the in- 
verse temperature (3. (Hereafter we take the unit K = 1.) 
The unit vector n represents the directions of staggered 



spin moments. The low-energy and the long-length scale 
behaviors of i?AFH are well described by NLerM as con- 
firmed by experiments 0] and quantum Monte Cairo 
simulations 0. 

There is a useful form for the NLerM, which is called 
CP 1 model 9] . Introducing complex fields Ca( x ) with 
= (c sw r,r) through n = £ CT , CT =T4 CowO with 
(a = x,y,z) representing Pauli matrices, we obtain 

S = (p s /2c sw )Jd 3 x [\d^( a \ 2 + (CM) 2 ] . Performing 
a Stratonovich-Hubbard transformation for the second 
term in the square brackets, we obtain the CP 1 model: 
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where g = c sw /2p° s and (a^(x)) = (i C ta {x)d^C ta {x)) 
at the saddle point. The field ev p is a U(l) gauge 
field associated with a local gauge symmetry, (, a (x) — > 
(a(x) exp(ix(x)). Thus, the NLerM is equivalent to the 
boson system with the U(l) gauge field interaction. 

The CP 1 model is also derived from the Schwinger 
boson mean field theory (SBMFT)[l(|. In the the- 
ory, the spin is represented in terms of boson op- 
erators Zj a and Zj a as Sj = ^J2 a ~' z ^" (Tr " J ' z 

with the constraint X) ct =t l z ]<j z j<? = 
SBMFT Hamiltonian is given by 

-<y/2) E(i,j) [ A * (^t z ji - z ^ z i\) + h - c - - 

X^jcr ^ ( z jo- 2 jo' — lj i w hh A and A being the mean field 
parameters. Here A = (zj, jZj \ — zn zg) describes the 
Schwinger boson pairing[lj,[lj,[lj|. After Fourier trans- 
forming, the Schwinger boson quasi-particle excitation 
spectrum is obtained by a Bogoliubov transformation^]: 
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uj k = J A 2 — J 2 |A| 2 (sin k x ± sinfcy) , where the plus 
sign is for k x k y > and the minus sign is for k x k y < 0. 
In the Neel ordering phase, A — 2J|A| and the Schwinger 
bosons are gapless at zone face centers (±7r/2, ±7t/2). 
The Neel ordering state is a consequence of Bose- 
Einstein condensation of the Schwinger bosons at those 
points ll4lll5l . One can derive the CP 1 model after some 
algebra 11]. The relation between the Schwinger boson 
fields and Ctr( x ) and (a(x) in Eq.@ is Zj a — Q a at one 
sublattice and Q a = zJ CT at the other sublattice. One 
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Schwinger boson carries spin 1/2 and the unit gauge 
charge. Therefore, pairs of the Schwinger bosons that 
describe the low-lying excitations carry integer spins and 
the gauge charge two. 

Both of NLcrM and SBMFT are, so to speak, bosonic 
description of the system. Meanwhile, one can construct 
a fermionic theory by representing the S = 1/2 spins 
in terms of fermions. Mean field analysis was carried out 
by Affleck and Marston[lj| and the 7r-flux phase was pro- 
posed. The state is characterized by a flux it penetrating 
each plaquette with alternating directions. The quasi- 
particle excitations are gapless at (±tt/2, ±tt/2). More 
advanced analysis based on an effective theory that in- 
cludes fluctuations around the mean field state through 
a U(l) gauge field shows that a mass term is induced in 
the spectrum dynamically. This phenomenon is associ- 
ated with the Neel ordering pL7L fia | . 

It has been shown experimentally that doped holes 
reside primarily on oxygen sites. Zhang and Rice 
suggested^^ that there is strong correlations of forming 
a singlet pair between oxygen p-orbital holes and cop- 
per d-orbital holes. The t-J model was proposed based 
on this picture. Although the single hole problem has 
been studied mostly by this model, the focus is mainly 
on frustration effects induced by hopping of the Zhang- 
Rice singlet [i^, EI- However, not so much attention has 
been paid on the spin configuration around the Zhang- 
Rice singlet. 

In this paper, it is argued that in the single hole doped 
system the doped hole induces a half-skyrmion which 
is schematically shown in Fig. ^ The half-skyrmion is 
a spin texture characterized by a half of a topological 
charge. 

We start with considering the spin configuration 
around a static Zhang- Rice singlet. The moving case 




FIG. 1: The half-skyrmion spin texture. Arrows indicate the 
staggered moment and the filled circle at the center represents 
the Zhang-Rice singlet formed site. 



shall be considered later. A Zhang-Rice singlet is com- 
posed of a d-orbital hole state and the four oxygen hole 
states around the copper ion. Constructing the Wan- 
nier functions, the operator for the Zhang-Rice singlet is 
given by 0| (<pj^dj[ — <j)j±djf), where is the sym- 
metric combination state of the four oxygen hole states 
around the copper ion at the site The d-orbital 

state of the Zhang-Rice singlet is given by the superposi- 
tion of the up-spin state and the down-spin state. To be 
specific, we assume that in the Neel ordered state before 
hole doping the ordered spin moment at the site j is in 
the direction of the positive z-axis. Then, the d-orbital 
state with up-spin does not affect the neighboring spins 
so much. 

In contrast, the down-spin state affects the neighbor- 
ing spins to change their directions. From the analysis 
of the NLcrM below, we show that a spin configuration, 
which is characterized by a non-zero topological charge, 
is created around the down-spin state. This spin con- 
figuration is called skyrmion[9j- For the static case, the 

energy is given by E = J d 2 r (Vn) 2 . The field equa- 
tion is V 2 n — (n • V 2 n) n = 0. Solutions are devided 
into sectors characterized by a topological charge : Q = 
i / (Previa ■ (<9 M n x d v n), where e xy = -e yx = 1 and 
txx = £yy = 0. The energy in each sector has the lower 
bound: E > A-Kp° s \Q\. The equality is satisfied if and only 
if d^n = ±e Mt , (n x d v n) 22] . The lowest enegy state is 
obtained by solving this equation with the boundary con- 
ditions: n — * +e z at infinity and n = — e z at r = Yj. In 
terms of a variable w = (n x + in y )/(l —n z ), the equa- 
tion turns out to be the Chaucy-Riemann equations |22(. 
Therefore, w is an analytic function of z — x + iy or 
z* = x — iy. The lowest energy state satisfying the 
boundary conditions is w s — [(x — Xj) + i(y — yj)]/X and 
w as = [(% — Xj) — i(y — Uj)]/X with A being a constant. 
In the vector n representation, these solutions are 




2X(x-x j ) 2X(y-y 3 ) jr-r^f-A 2 ^ 



(3) 

where the plus sign is for uj s and the minus sign is for 

Was- 

Now we consider the superposition of the uniform state 
n = +e 2 for the up-spin state and the skyrmion spin con- 
figuration © created by the down-spin state. The NLcrM 
analysis suggests that a spin configuration with non-zero 
topological charge Q, where < \Q\ < 1, is formed with 
ni z > for any site £. The value of the topological charge 
Q is determined through the flux representation in terms 
of the CP 1 gauge field a u : Q = ± J d 2 r (V x a) . The 
spin configuration with Q corresponds to the flux 27r<3. 
Due to Bose-Einstein condensation of the bosons C<r( x )> 
the value of the flux 27rQ is not arbitrary. Since low- 
lying excitations are pairs of these bosons, the flux quan- 
tum is 7T analogous to conventional BCS superconduc- 
tors. From the constraint on Q, i.e., < \Q\ < 1, we 
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conclude that 27r|Q| = tt, that is, \Q\ = 1/2. There- 
fore, the spin configuration is the half-skyrmion, which 
is equivalent to the skyrmion with the core, and its en- 
ergy is E = 2irp s . This conclusion can be reached if we 
assume that the staggered magnetization vanishes at the 
perimeter of the core. The analysis of the NLctM gives 
the half-skyrmion solution as was shown by Saxena and 
Dandoloff [23j • They considered the two-dimensional fer- 
romagnetic Heisenberg model in the context of the quan- 
tum Hall systems. For the static case, the field equation 
is the same. So we can apply their analysis to the anti- 
ferromagnetic case as well. 

The moving half-skyrmion is obtained from the static 
half-skyrmion by applying a Lorentz transformation be- 
cause the action JU and J5J are Lorentz invariant 
with c sw the speed of "light." From the calcula- 
tion of the energy-momentum tensors, we find e° = 

^c 2 w (/c 2 + A: 2 ) + E$. For the S = 1/2 antiferromag- 

net, the bare spin- wave velocity is c sw = \[2J and the 
bare spin stiffness is p® = J/4. Renormalization ef- 
fect leads to c sw — * \J~2Z C J and p a s — > Z p J/4(= p a )\$- 
Thus, the renormalized dispersion is given by e& = 
J^2Zl(kl + fc 2 ) + (ttZ p /2) 2 . Since the doped hole is 

fermion, this relativistic dynamics of the half-skyrmion 
is described by a Dirac fermion Lagrangian density: 

£ = X^( 7 ^ + mc sw)Vv, (4) 

with tp = V-^7o and toc 2 w = irZ p J/2. The 7 matrices sat- 
isfy {7^,7^} = 2S plJ in the Euclidean space time. Here 
a is an index for the sign of the topological charge. The 
number of components of the Dirac fermion field is four 
because there are a positive energy state and negative en- 
ergy state and the origin of the Dirac fermion dispersion 
is at either ki = [tt/2,tt/2) or k 2 = (— k/2, tt/2). The 
latter is implied from the fact that the Schwinger bosons 
are gapless at those points in the symmetry broken Neel 
ordering phase. Note that the dynamics is associated 
with the properties of the parent compound, namely, the 
Mott insulator. In this sense, the fermion ip is different 
from quasi-particles in the conventional Fermi liquid. 

The lattice version of the action is derived by discretiz- 
ing Eq.EJ. As is well-known in the lattice gauge field 
theory [24j, there is a fermion-doubling problem when one 
introduces a lattice for a Dirac fermion. This subtlity 
does not make any difficulty in our case because of the an- 
tiferromagnetic correlations that naturally introduce two 
species of Dirac fermions with half of the Brillouin zone 
(the magnetic Brillouin zone). Discretizing Eq.QJ and 
then performing Fourier transformation, we obtain C = 
„ -j ( mc 2 w + d T cos k x + i cos k y \ 

l^a Vka y _ cog kx + i CQS ky mc 2 w _ Q r j Vka ■ 

Here we have used the fact that the origin of the Dirac 
fermion dispersion is at ki2- The dispersion relation is 



given by 



e fe = ± \ c sw( cos2 k x + cos 2 k y ) + (toc sw ) 2 . (5) 

The dispersion curve ej~ — toc 2 w is qualitatively in good 
agreement with the experimentally observed dispersion 
by Wells et al^, . In our theory, Z c and Z p are deter- 
mined by the parameters of the undoped system if one 
neglects spin wave effects on the self-energy of the half- 
skyrmion. For the S — 1/2 antiferromagnetic Heisenberg 
model, the values of Z c and Z p are estimated as Z c = 1.17 
and Z p — 0.72 from quantum Monte Cairo simulations 
and a series expansion analysis |2^|. Substituting these 
values into Eq. we find that the band width along 
(0, 0) to (tt, tt) is ~ 1.5 J and the Dirac fermion mass is 
to ~ 1.13J. Experimentally, the band width is estimated 
to 2.2J 3]. This 30% discrepancy would be associated 
with the deviation of the real system from the NLctM 
description. There is also self-energy corrections on the 
mass to due to spin waves. We have estimated it to 
to/ J = 1.13/ (1 — 0.12e^) where is a coupling parame- 
ter between the half-skyrmion and the antiferromagnetic 
spin waves. The coupling parameter turns out to be 
a gauge charge in the dual theory below. 

The effective action for the half-skyrmion can be de- 
rived by duality mapping [2?|. In the CP 1 theory, the 
half-skyrmion solution residing at the origin has the fol- 
lowing form z a (x) = Xu "+ 1 J^™-p( ±ie ) ; w here u* • v = 0, 
and |u| = |v| = 1. (For the solution J3J with set- 
ting rj = 0, these parameters are 77 = U[ = 1 and 
M| = = 0.) If we take \v^\ — |wj ^ 0, then the stag- 
gered moments lie in the plane with the same direction at 

infinity. In this boundary condition, z a (x) ~ v a exp(±z#) 

1/2 

at r ^> A. We write z„(x) — pj exp(z0(x)), where 
4>(x) — 4>o(x) + 4> v (x). Here <f>o(x) is a non-singular func- 
tion of x and 4> v (x) describes the vortex (half-skyrmion) 
field: 4> v (x) ~ q v tan _1 (y — y v )/(x — x v ), where q v is the 
sign of the topological charge and (x v ,y v ) denotes the po- 
sition of the vortex. (Note that x v and y v are functions 
of the imaginary time.) In applying the duality mapping, 
a dual gauge field is introduced. The coupling between 
the gauge field and the half-skyrmion has the form of the 
minimal coupling. Thus, we obtain the following QED3 
Lagrangian density in the continuum, 

C = ^2 [-f p [dp - iq a A^) + to] {dpAv - d u A p ) 2 . 

(6) 

The gauge field originally comes from the spin cur- 
rent, which is associated with the antiferromagnetic spin 
wave modes. The gauge charge ca is a parameter of the 
theory. Estimation of is beyond the analysis of the 
effective theories. For analysis of the half-skyrmion self- 
energy, we formulate the theory on the lattice according 
to a standard procedure [24|. The second order term with 
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respect to the one "photon" vertex part is evaluated nu- 
merically and the first order term with respect to the two 
"photon" vertex part is computed analytically. While 
the former has negligible effect, the latter leads to — 
-(^ei/16) 7, sin M 

The mass renormalization due to this term is 
m -> m/(l - 0.12e^). 

A similar action of Eq. © is obtained in the fermionic 
representation of the 5=1/2 antiferromagnetic Heisen- 
berg model 0, 0] based on the 7r-flux phase^ with 
dynamically induced mass m, through the interaction 
with a gauge field. (The condition of this dynamical 
mass generation 28J is that the number of Dirac fermion 
species is less than 32/?r 2 ~ 3.2. This is the case for the 
QED3 theory of the ant if erromagnet 1 1 & 12 91 . ) From exact 
diagonalization studies it was shown [3(| that the bond 
spin currents which characterize the 7r-fl ux p hase is re- 
produced by a skyrmion-like spin texture |31l. From the 
view point of the spin-charge separation 32], Baskaran 
suggested!^ that a half-skyrmion can be seen as de- 
confined spinons based on an analysis of an 71-skyrmion 
solution. 

In summary, we have argued that, within the effective 
theory approaches to the S = 1/2 Heisenberg antiferro- 
magnet, the doped hole induces the half-skyrmion spin 
texture through the formation of the Zhang-Rice singlet 
in the single hole doped system. This picture is consis- 
tent with rapid suppression of the Neel ordering by the 
doped holes because the half-skyrmion behaves as a vor- 
tex in the condensate. It would be interesting to extend 
the picture to the slightly doped regime. 
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